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When a liquid is cooled below its melting temperature it usually crystallizes. 
However, if the quenching rate is fast enough, it is possible that the system re- 
mains in a disordered state, progressively losing its fluidity upon further cool- 
ing. When the time needed for the rearrangement of the local atomic structure 
reaches approximately 100 seconds, the system becomes "solid" for any practi- 
cal purpose, and this defines the glass transition temperature Tg. Approaching 
this transition from the liquid side, different systems show qualitatively dif- 
ferent temperature dependencies of the viscosity, and, accordingly, they have 
been classified introducing the concept of "fragility". We report experimental 
observations that relate the microscopic properties of the glassy phase to the 
fragility. We find that the vibrational properties of the glass well below Tg are 
correlated with the fragility value. Consequently, we extend the fragility con- 
cept to the glassy state and indicate how to determine the fragility uniquely 
from glass properties well below Tg. 

When a liquid is cooled, the loss of kinetic energy leads to an ordering of the molecules 
which then crystallize at the melting temperature T^. However, if cooled fast enough through 
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Tm, some materials (glass forming materials) are capable to sustain a metastable liquid state and, 
upon further cooling, to freeze into a disordered glassy state (|7]|2]|5llll). The law that describes 
the change of the viscosity with the temperature approaching the glass transition temperature, 
Tg, is highly material specific, and has led to the classification of the glass formers materials 
according to the concept of "fragility" (l5|0). The kinetic fragility, m, is directly related to the 
slowing down of the dynamics: it is defined in terms of the shear viscosity 77 as: 

d\og{ri) 

m = lim 

T^T, d{TjT) 

Therefore, m is an index of how fast the viscosity increases approaching the structural ar- 
rest at the glass transition temperature Tg, the temperature where the structural relaxation time 
Tq ~ 100 s. At this latter temperature, through the Maxwell relation rj = CqoT, corresponds a 
viscosity r] ~ 10^^ poise (more likely i] ~ 10^^ poise for molecular glasses), while 10^"^ poise 
is the "infinite" temperature limit in basically any materials. Consequently, fragility values typ- 
ically range between m = 17 for "strong" systems -those that show an Arrhenius behaviour- 
and m ~ 150 for "fragile" systems, where the high cooperativity of the diffusive dynamics 
induces a high (and T-dependent) apparent activation energy. One interest in this classification 
lies in the attempt to relate the temperature behavior of a macroscopic transport property close 
to Tg to the microscopic interactions driving the dynamics of the system. It has been found, for 
example, that the value of the fragility is empirically related to the kind of interaction potential 
among the particles constituting the system. Prototypical examples of fragile liquids are those 
composed by units interacting via isotropic bonds, such as Van der Waals-like molecular liquids. 
The strong glass-forming liquids, on the other hand, are those characterized by strong covalent 
directional bonds that form space filling networks. 0-terphenyl (m=80) and 5*^02 (m=20) are 
characteristic examples of a fragile and a strong liquid, respectively. Hydrogen bonded systems, 
as glycerol (m=50), are often called "intermediate" between strong and fragile liquids. The ki- 



netic fragility has been found to be correlated to other properties of the glass-forming liquids, 
such as: i) the slope at Tg of the temperature dependence of the configurational entropy (often 
referred to as thermodynamic fragility) ((Tj), or -for classes of systems sharing similar glass tran- 
sition temperatures- the specific heat jump at Tg (013); ii) the ratio between the maximum and 
the minimum of the boson peak, i.e. of the bump observed at the THz in the Raman and neu- 
tron scattering spectra of glass-forming materials (17^ . but this finding is controversial (1771 : iii) 
the degree of stretching in the non-exponential decay of the correlation functions in the liquid 
close to Tg {\12ii . iv) the statistics of the minima in a potential energy landscape-based descrip- 
tion (|7^|7?|) of the diffusion process in supercooled liquids (|75l|73|73), and, more recently, v) 
the temperature behaviour of the shear elastic modulus in the supercooled liquid (I78t . In all 
these studies the fragility has been always related to (or defined through) macroscopic proper- 
ties characterizing the liquid side of the glass -transition. While there are attempts to relate the 
fragility to the anharmonicity of the "hot" glass (179b . no connection has been found up to now 
between the value of m and the physical properties of the low temperature glassy phase. 

We show that, starting from a determination of the non ergodicity factor in the low temper- 
ature glass, it is possible to identify a parameter that controls how fast the non ergodicity factor 
decreases on increasing the temperature, and that turns out to be proportional to the fragility 
m. Through this, we establish a way to determine the fragility of a system in the glassy phase 
well below Tg, independent of the way the viscosity changes with decreasing temperature from 
the liquid side. By exploiting the harmonic approximation of the low temperature dynamics, 
it is found that this parameter only depends on the characteristics of the static disorder, which, 
in turn, controls the vibrational eigenmodes of the glass. This result demonstrates the exis- 
tence of a deep link between the diffusive inter-basins dynamics and the vibrational intra-basin 
dynamics. 

Recent extensive inelastic x-ray scattering (IXS) measurements of the dynamic structure 
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factor have allowed to constitute a sizeable library of high frequency (THz) dynamical prop- 
erties of glasses. Of interest here, the IXS measurements allow for the determination of the 
non-ergodicity factor, f(Q,T) with a reliability that was not achievable via light or neutron 
scattering S2U\i . The non ergodicity factor is the long time limit of the density correlator, $Q(t), 
i. e. the density-density correlation function, F{Q, t), normalized to the static structure factor, 
S{Q): $Q(t) = F{Q, t)/S{Q). The quantity l-f(Q, T) represents the amount of decorrelation 
introduced by the vibrational dynamics, and it depends on both the (T-dependent) amplitude 
of the vibrations and the degree of disorder of the glassy structure. In a low temperature glass, 
F(Q,t), apart from the Debye- Waller factor exp {—W{Q)), can be expressed as the sum of 
a constant term Sis{Q), which represent the static structure factor of the atomic equilibrium 
positions (Inherent Structure), plus a time-dependent one, FimiiQ, t), which is the contribution 
of the atomic vibration around such equilibrium positions, a quantity that vanishes in the long 
time limit: 

F(Q, T) = e-^(«) [SisiQ) + F^nei{Q. t)] (1) 

Therefore: 

where we have defined SineiiQ) = FineiiQ, t = 0). S{Q, cu) is the Fourier transform of F{Q, t) 
and is the quantity directly accessible in scattering experiments. From Eq.([l]) it can be expressed 
as 

SiQ, u) = e-^(«) [SisiQ)5icu) + S,UQ, cu)] (3) 
From an experimental point of view and according to Eqs Q and Q, the non ergodicity 
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factor is derived from the ratio of the elastic to the inelastic scattered intensity, obtained from 
inelastic scattering measurements of the dynamic structure factor S{Q, u) (|27]l . A sense of the 
T dependence of f{Q, T) can be obtained from Fig. 1. Here, as an example, we report the IXS 
spectra at fixed exchanged wavevector {Q = 2 nm^^) and at different temperatures in glycerol. 
The inelastic (dashed lines) and elastic (dotted lines) contributions to the scattering intensity 
are shown and one can appreciate in the raw data the change of relative intensity as a function 
of T. As far as the Q dependence is concerned, f{Q,T) follows in phase the oscillations 
of the static structure factor and is almost Q-independent in the Q ^ region where S{Q) 
is almost constant (E^ . see the inset of Fig. 2. We focus on this small-Q region. From the 
integrated intensities of the elastic and inelastic contributions, obtained by a fitting procedure 
(see Eq. (SI) and (S2) of the supporting on line materials), the T dependence of f{Q,T) is 
obtained. The values of f{Q, T)~^ (which is expected to be linear in T, see below) are reported 
in Fig. 2 (triangles). Also reported in the same figure is the T dependence of /(Q, T) for two 
other archetypical glasses: silica and o-terphenyl (oTP). 

To better understand the temperature dependence of f{Q, T) in the T ^ limit, we invoke 
the harmonic approximation for the vibrational dynamics. This allows one to express f{Q,T) 
in terms of the vibrational properties of the systems, i. e. the eigenvalues (ujp) and eigenvectors 
(Cp) of the the potential energy Hessian evaluated at the inherent structure. Using the harmonic 
approximation for SineiiQ, tu), it is straightforward (I23II24I) to show that Eq.© reduces to: 

-1 

(4) 

Here M is the molecular mass, Kb the Boltzmann constant, and i is summed over the 
particles and p over the 3A^ normal modes. In order to pinpoint the T dependence of the non 
ergodicity factor in the low Q limit it is convenient to rewrite Eq. Q as 
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f{Q-^0,T) = —^. (5) 

1 + aTfr 

We thus define the dimensionless quantity a, which encompasses all the microscopic details 
of the system, as the eigenvalues and eigenvectors of the normal modes. These are quantities 
that, in turn, depend on the interaction potential and on the disordered structure. This equation 
provides a formal way to extract the system-dependent parameter a from the T dependence of 
f{Q, T), derived from the IXS data. This has motivated us to revisit the large amount of IXS 
data available for glasses at low T where the harmonic approximation, and therefore Eq.(|51), 
is expected to be valid. As it can be seen in Fig. 2, the observed T dependence of f{Q, T) is 
fully consistent with the functional form predicted by Eq.(|5l), and this allows us to determine 
a by a least square minimization procedure. The derived values for a (e. g. a=0.l9 for silica, 
a=0.32 for glycerol and a=0.58 for oTP) clearly indicate a trend: the more fragile a liquid, the 
greater the slope of f{Q, T) at T=0, i. e. the faster the decorrelation of the density fluctuations 
on increasing T. The fitting procedure has been applied to the whole set of available glasses, 
and the obtained values of a are reported in Tab.I and Fig. 3 as a function of the independently 
known fragility parameter m. From this figure clearly emerges the existence of a strong correla- 
tion between m and a: the higher the fragility, the higher the value of a, i. e. the faster is the T 
dependence of the f{Q,T) parameter. The existence of a strong correlation between a and m is 
further emphasized by the empirical observation that the two quantities are not only correlated 
but (within the statistical accuracy) actually proportional to each other, according to the relation 
m = (135 ± 10)a + (4 ± 5). On passing, we note that the two points that lie definitely below 
the dotted line (Selenium and Salol) are the ones for which fragility determined at Tg doesn't 
agree well with the fragilities determined at higher temperatures (171 E5b. 

The observed correlation is conceptually surprising. It indicates the existence of a link be- 
tween the curvatures of the potential energy function at its minima (more specifically those 
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visited in the glassy phase) and the other properties of the potential energy function (energy dis- 
tribution of the minima, minimum-to-minimum barrier height, distribution of the saddle order 
and energies,...) controlling the diffusive processes in supercooled liquids. 

We further examine how a emerges from the collective density-density correlation function 
plateau. From Eqs. ^ and dSj), the microscopic expression for a is found to be: 



One may in principle derive a similar parameter ag from the temperature dependence of the 
self correlator plateau. In this case is related to the familiar mean square displacement. 
The harmonic approximation for the Debye Waller factor /^(Q, T) = exp (—W{Q, T)) would 
lead to an equation formally identical to Eq.©, but with a replaced by a^: 



Therefore a and differently weight the low frequency modes. Specifically, in the small Q 
limit, a is more sensitive to the low energy modes than as that, independently on Q, reflects the 
whole density of states. It would be interesting to check whether or not these two quantities are 
correlated with each other. 

In conclusion, we report evidence for the correlation between the fragility of a glass-forming 
liquid and the temperature dependence of its non-ergodicity factor as determined by the vibra- 
tional dynamics at very low temperatures. The fragility is an index of how the viscosity in- 
creases upon supercooling. The non ergodicity factor, in the low temperature limit, is related to 
the vibrational properties of the harmonic glassy dynamics (see Eq.©), i.e. to the curvature of 
the energy minima. Therefore, from our finding it emerges that the properties of the potential 
energy landscape around the deep minima are related to those properties that control the diffu- 
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sion of the system through different basins. This unexpected relation represents a further aspect 
that requires to be clarified in the physics of the glass-transition. 
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Table 1: Temperature steepness of the viscosity at Tg (fragility) and of the non ergodicity factor 
at T ^ (a) for several materials. 

from Ref. a from Ref. (27). 
^Silica infrasil grade sample, m from Ref. (17^ . a from Ref. (EH) and references therein. 

from Ref. (|i2|), a from Ref. (29). 
'^Poly butadiene: a is determined from an experiment performed on a (1.2PBD)o(1.4PBD)ioo 
reported in Ref. (0^. The fragility is not available for such concentration, we estimated m 
extrapolating data from Refs. i31\U2i . 

^ Normal-butyl-benzene: m from Ref. (0^ . a from the experiment in Ref. 

from Ref. (f72l) . a from unpublished data. 
^Meta-toluidine: m from Ref. (34), a from unpublished data. 

Orto-terphenyl: m from Ref. (72), a from Ref. SJSl . 
* Meta-tricresyl-phosphate: m from Ref. (122), a from unpublished data, 
from Ref. dJ^ . a from unpublished data. 
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Figure 1: Selected example of the inalastic x-ray scattering energy spectra of glycerol (open 
circles with error bars) taken at Q = 2 nm^^ at the indicated temperatures. The solid line is the 
line of best fit according to Eqs ^ and (S2), while the dashed and dotted lines are the elastic 
and inelastic contributions, respectively (see the supporting on line material for further details). 
The values of f{Q, T) are obtained by the ratio of the integrated intensities of the elastic and 
inelastic contribution, from Eq.©. ^3 
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Figure 2: Values of f{Q*, T)-^ (Eq.©) for Q* ^ 2 nm-^ reported in a T/Tg scale for three 
representative materials (full symbols with error bars): silica (Tg=1450 K), glycerol {Tg=\90 
K) and oTP (Tg=241 K). The full line is the best fit of the experimental data to Eqs These 
fits have been used to derive the values of a reported in Fig. 3. In the inset we show the 
Q-dependence of /(Q, T*) for Silica at T* = 1050 K 
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Figure 3: Correlation plot of the kinetic fragility and the a parameter of the non-ergodicity 
factor (see Eq. (jSj). The dotted line is obtained by a fit of the data to a linear equation. It 
corresponds to m = 135a; and the regression coefficient is r=0.994. 
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